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In this paper the polaron problem for the Holstein model is studied in the weak coupling limit. We 
use second order perturbation theory to construct renormalized electron and phonons. Eigenstates 
of the Hamiltonian are labelled and the excitation spectrum is constructed. 

pacs numbers: 71.38.+i, 63.20.Kr 



> 



OC 



o 
o 



00 ■ 

Q\ • I. INTRODUCTION 

on : 

The Holstein model which describes electrons coupled to some local molecular deformations is a simple model 
Oh. to study polaron properties (QJ2). These recent years, accurate numerical investigations of this model on finite size 
systems have given new insights closely linked to some current problems in condensed matter physics |2j. A 

polaron corresponds to an electron renormalized by the interactions with the lattice deformations or phonons || . From 
a historical point of view |g|-|Tl|l the approach to the polaron problem consists in computing the spectrum of the whole 
system containing a single electron coupled to some optical phonons. Writing k the momentum of the electron and 
Q the momentum of the phonons, the problem is to compute the lowest energy of the system for a total momentum 
K = k + Q, since the total momentum operator commute with the Hamiltonian. This mimimum eigenenergy of the 
Hamiltonian for a given total momentum K gives the polaron relation dispersion, as far as we consider only the low 
\Q lying excitations of the system -that is as far as no phonon excitations are involved. 

7— I , In this paper, we consider a different approach to the problem and try to decompose the excitations of the system 
0^ ' in term of polaron excitations and renormalized phonons excitations. Due to the difficulty of the problem, we address 
only the weak coupling limit of the problem and we use mainly second order perturbation theory in the electron 
phonon coupling. Within this framework, one can block diagonalize the Hamiltonian. Then it is possible to label 
the eigenstates of the Hamiltonian in term of renormalized electron and phonons. In other words, an eigenstate of 
the Hamiltonian contains one renormalized electron (a polaron) of momentum k and some renormalized phonons of 
momentum Q; the total momentum of the state being K = k + Q. 
""j" 1 , The Hamiltonian for the Holstein model, in Wannier representation, is given by 
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g ; H = -tYc], sCj + woYbfa - gwoYctefa + bl), (1) 
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where c\J and Cj are electron creation and annihilation operators, bj and bj are phonon creation and annihilation 
operators, t is the hopping integral, loq is the optical phonon frequency and g is the dimensionless coupling constant. In 
the strong coupling limit, g 1, the Hamiltonian is approximatively diagonalized by the Lang-Firsov transformation 
|^2|,p|. The transformed Hamiltonian is given by H — e~ s He s , with S = J2j 9 c j c j{bj~bj)- It is only when the hopping 
term is zero that the transformed Hamiltonian is diagonal. In this case one obtains some new fermionic operators, 
which describe small polarons, or localised polarons. The small polaron creation operator is given by pt = e s cje _s , 

or pj — UjCj, with Uj = exp g(b| — bj) while the renormalized phonons operators are given by Bj = bj — g. The 
action of this small polaron creation operator on the vacuum is thus to create an electron and a coherent state for the 
phonons where the annihilation operator takes the expectation value g. For non zero values of the hopping integral, 
these small polarons become some quasiparticles more or less well defined depending on the value of the coupling and 
on the phonon frequency |^,|l3|]. We notice that for g = 0, the small polaron correspond to the bare electron, since 
there is no coupling and thus no renormalization. We shall use a similar unitary transformation formalism but in we 
weak coupling limit, g — > 0, and work in momentum space where the Hamiltonian reads 

H = + ^oY b % - T^EcU^ + fo^) (2) 

fc q ^ k,q 

where M is the number of sites of the lattice, and indices of various operators correspond to Bloch momenta. In the 
next section we consider our approach for a two site system, then we will generalize to a general system. 
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II. THE TWO SITE SYSTEM 



We start our perturbation theory for the two site system. In this case, the momentum takes two values k — 0, tt 
and the bare electron relation dispersion is given by Sq = —t and = +t. The q = mode of the phonons is 
coupled to the total number of electron which is one in this problem. We drop out this mode after performing the 
shift corresponding to the displaced oscillator 0. The reduced Hamiltonian is then given by 



He 



(cjc + 4 C *0 - t (cjc - C^Ctt) + w tfb 



guo ( t 

V2 " 



(c c 7r + ctc )(fc + fe t ) 



(3) 



We have dropped the index tt for the boson operators. The basis states are given by \k\ n) — c\(tf ) n / VnJ.\0) and |0) is 



the vacuum. For g — the corresponding eigenenergies are given by e n {k) 
perturbation theory, we obtain the eigenstates and eigenenergies, 



nw + £fc — g Wq/2. Using second order 



E n (k) = e n (k) + 
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2ek + cop 



+ (n + l) 
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(4) 



\k;n) 



gup 
V2 



2t — ujq 2t + wo 
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2t + wo 2t — wo 
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\k\n) 



(5) 



We show now that the eigenstate \k;n) describes a state with one polaron of momentum k and n renormalized phonons. 



The states \k,n) are eigenstates of the Hamiltonian Hr, thus using Hr 
to 0(g 2 ), the S- matrix of the unitary transformation is given by 



'ii^e , we obtain \k;n) 



gup 
V2 



6t 



2t — ujq 2t + wq 
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2t + Wq 2i — Wq 



|fe;n). Up 



(G) 



We note that if we set t = in this expression, we obtain the right transformation which diagonalize the reduced 
Hamiltonian Hr. We can now construct the corresponding polaron operator pk through the transformation p k — 



e s c k e 



~ s . We obtain explicitly 



Po = c 



g^o 
V2 



2t + Wq 2t — Wq 



(7) 



gujo 
V2 



2t — Wq 2£ + Wq 



(8) 



These operators are fermionic (they satisfy the anticommutation relations) , with the property p\pk \k';n) — Sk.k'\k',n). 
This means that the eigenstate \k; n) contains one polaron of momentum k. As concern the phonons, we can proceed 
the same way, introducing a new operator B = e s be~ s which describes a renormalized boson or phonon. Explicitly, 
one obtains 



V2 



6t 



2t — Wq 2t — Wq 
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(9) 



Again, this operator is diagonal with B^B\k;n) — n\k;n). We therefore conclude that the state \k;n) contains one 
polaron of momentum k and n renormalized phonons. The total momentum of the state is K = k if n is even and 
K = k + n if n is odd. Then we write the eigenenergies as E n (k) = e* ± t* . One obtains 



tiwq 



<7 2 w U 2 - 2uj 2 



At 2 



(10) 



t 1 



wn 
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For n phonons or renormalized phonons in the system t* corresponds to the fermionic (polaronic) excitation energy. 
In the polaron problem, we restrict ourselves to the case n = 0, that is we create a polaron from the vacuum, 
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|fc; 0) = p k \0). We are left with e* = Eq which is the polaron binding energy and t* = t$ which represents the hopping 
integral of the polaron. If we take g = 0.4 and uiq = 10/llt we obtains e* = — 0.0538i and t* = 1.042i while exact 
diagonalization || gives s* = — 0.0576£ and t* = 1.038i. We notice the enhancement of the bandwidth for this two 
site system. 

In term of the renormalized operators, the Hamiltonian can now be written in a block diagonal form, 

H = J2 E o(k)p{ P k + ^ -i^oWl^U 
k ' k 



+ y Y} E ^ - 2E ^ + Eo(k)]Bt&BBplp k + ... (12) 

k 

The original electron phonon interaction which contained off diagonal transitions between subspaces with the same 
number of particle (electron or phonon) has been eliminated: there is no more off diagonal interactions between the 
polaron and the renormalized phonons. 



III. GENERAL CASE 



Next, we consider finite size systems for which numerical results have been obtained 



the configuration of phonons {ft} such that b*b q \ft 



The basis vector are then \k:n) 



Let In) be the state with 



The eigenstates 



ct\n) 



are given by \k;ft) = e s \k; n) 



Now we use the fact that the state \ft) contains no electron, so that e s \n) — \ft). Using 
e s cke~ s , we obtain the desired relation, p\\ft) — \k;ft). In other words, 



the definition of the polaron operator pk 

the operator p\ acting on the state \n) with the configuration {ft} of bare phonons, gives the eigenstate \k\ n) which 
contains one polaron of momentum k and the configuration {ft} of renormalized phonons. Defining B q = e s b q e~ s , we 
therefore obtain B^ q B q \k\ ft) — n q \k; ft). This is the straightforward generalization of the results obtained for the two 
site system. 

Using second order perturbation theory the polaron dispersion is given by 



E(k) = e k 



2 2 

M 
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(13) 



This is a well known result obtained by many approaches (^J^). This dispersion relation corresponds to a state with 
one polaron of momentum k and zero renormalized phonons. The next excited states correspond, for example, to 
states with one polaron of momentum k and one renormalized phonon of momentum q. The total momentum of these 
states is therefore K = k + q. The energy of such a state is given by 



E(k;q) - E(k) 
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_£fc — £k+q + ^>0 £k — £k-q — 

As concern the corresponding unitary transformation, it has been discussed by Frohlich |14| and is given by 



S 



5^0 \ - 



hi 



M ^ 
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£fc + ^0 — e k+q £k ~ £k+q — ^0 



"l+q C 



(14) 



(15) 



We now discuss the excitation spectrum of the system in the weak coupling limit and consider a one dimensional 
system for simplicity with the bare electron dispersion Ek = — 2tcos(fc). 

Our basic asumptions are a) that the dispersion of the polaron E(k) is slightly renormalized from the bare electron 
dispersion relation Ek and b) the system is large enough so that the phonon frequencies are slightly renormalized, so 
that E(k;q) ~ E(k) + luq. The ground state corresponds to the state with a single polaron of momentum k = 0, 
with energy E(0). The state of momentum K with the lowest energy can be either the state with one polaron of 
momentum k = K and no renormalized phonons or either the state with one polaron of momentum k = and one 
renormalized phonon of momentum q = K such that K = k + q. The corresponding energies are respectively E{K) 
and Eifi) + cjq. Within our asumptions on E(k), there exists a threshold K* such that E(0) + luq = E(K*). For 
K < K* the lowest energy is E(K) and the state contains only one polaron of momentum k — K. For K > K* 
the lowest energy is E(0) + cjq and the state contains one polaron of zero momentum and one renormalized phonon 
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of momentum q — K. As far as the others excitations are concerned, their construction is straightforward. If we 
compare this spectrum with the results of Ref. 0j^], where this spectrum was computed, we find a nice agreement 
in the case of a weak coupling. The flat part in the spectrum is just the energy E(0) + u>q. We further notice that 
for K < K* , these numerical results show that the polaron dispersion E{k) is slightly renormalized from the bare 
electron dispersion 

IV. CONCLUSION 

In conclusion, we have shown that for any coupling, one can define a renormalized electron of momentum k and 
some renormalized bosons of momentum q, such that the spectrum is made up of these renormalized excitations. The 
renormalized operators can be obtained via an unitary transformation. This was worked out explicitly in the weak 
coupling limit using perturbation theory. Recent finite size studies support our approach. 

In the crossover regime, where level crossings arise in the spectrum of the Hamiltonian , it is not possible to define 
the unitary transformation pertubatively and it is not clear how to make some quantitative predictions. 
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FIG 1. The dispersion of the polaron for a one dimensional lattice with uja = 0.8 and g 2 u>o = 0.1. Energies are in unit of t, 
solid line is the bare electron dispersion, dotted line is the slight polaron dispersion and dashed line is the lowest excitation as 
discussed in the text. 
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